REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  V A  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 


reprint 


4.  TITLE  AND  SUBTITLE 

Transfer  of  rovibrational  energies  in  hydrogen  plasma-carbon  surface 
interactions 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 

W91  INF-05-1-0265 

5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

P.S.  Krstic,  E.M.  Hollmann,  C.O.  Reinhold,  S.J.  Stuart,  R.P.  Doerner,  D. 
Nishijima,  A. Yu.  Pigarov 


5d.  PROJECT  NUMBER 

W91  INF-05-1-0265 

5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Thompson  &  Sewell  Research  Groups,  Department  of  Chemistry 
University  of  Missouri-Columbia 
Columbia,  MO  65211 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.  S.  Army  Research  Office 
P.O.Box  12211 

Research  Triangle  Park,  NC  27709-2211 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  federal  purpose  rights. 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

48101-EG-MUR 


14.  ABSTRACT 

We  present  state-of-the-art  molecular  dynamics  (MD)  simulations  of  high-density  plasma-bombardment  of  carbon,  mimicking 
plasma-wall  interactions  at  the  fusion  reactor  first  wall.  Bare  and  hydrogenated  amorphous  carbon  surfaces  with  temperatures  in  the 
range  of  300-800  K  are  bombarded  by  a  distribution  of  neutral  hydrogen  molecules  representing  well-defined  center-of-mass  and 
rovibrational  temperatures.  The  MD  simulations  are  benchmarked  against  experiments  in  which  a  heated  carbon  surface  is 
irradiated  with  hydrogen  molecules  from  a  plasma  source.  Comparisons  between  simulations  and  experiment  are  presented  for  the 
rovibrational  distributions  upon  reflection  and  the  rotational  and  vibrational  accommodation  coefficients. 


15.  SUBJECT  TERMS 


chemical  sputtering,  amorphous  carbon,  accommodation  coefficients,  plasma 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

U 

U 

U 

u 

4 

Dr.  Donald  L.  Thompson,  PI 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

573-882-0051 


Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


INSTRUCTIONS  FOR  COMPLETING  SF  298 


1.  REPORT  DATE.  Full  publication  date,  including 
day,  month,  if  available.  Must  cite  at  least  the  year 
and  be  Year  2000  compliant,  e.g.  30-06-1998; 
xx-06-1998;  xx-xx-1998. 

2.  REPORT  TYPE.  State  the  type  of  report,  such  as 
final,  technical,  interim,  memorandum,  master's 
thesis,  progress,  quarterly,  research,  special,  group 
study,  etc. 

3.  DATES  COVERED.  Indicate  the  time  during 
which  the  work  was  performed  and  the  report  was 
written,  e.g.,  Jun  1997  -  Jun  1998;  1-10  Jun  1996; 
May  -  Nov  1998;  Nov  1998. 

4.  TITLE.  Enter  title  and  subtitle  with  volume 
number  and  part  number,  if  applicable.  On  classified 
documents,  enter  the  title  classification  in 
parentheses. 

5a.  CONTRACT  NUMBER.  Enter  all  contract 
numbers  as  they  appear  in  the  report,  e.g. 

F3361  5-86-C-51  69. 

5b.  GRANT  NUMBER.  Enter  all  grant  numbers  as 
they  appear  in  the  report,  e.g.  AFOSR-82-1  234. 

5c.  PROGRAM  ELEMENT  NUMBER.  Enter  all 
program  element  numbers  as  they  appear  in  the 
report,  e.g.  61 1  01  A. 

5d.  PROJECT  NUMBER.  Enter  all  project  numbers 
as  they  appear  in  the  report,  e.g.  1  F665702D1 257; 

I  LI  R . 

5e.  TASK  NUMBER.  Enter  all  task  numbers  as  they 
appear  in  the  report,  e.g.  05;  RF0330201 ;  T41 12. 

5f.  WORK  UNIT  NUMBER.  Enter  all  work  unit 
numbers  as  they  appear  in  the  report,  e.g.  001; 
AFAPL304801 05. 

6.  AUTHOR(S).  Enter  name(s)  of  person(s) 
responsible  for  writing  the  report,  performing  the 
research,  or  credited  with  the  content  of  the  report. 
The  form  of  entry  is  the  last  name,  first  name,  middle 
initial,  and  additional  qualifiers  separated  by  commas, 
e.g.  Smith,  Richard,  J,  Jr. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND 
ADDRESS(ES).  Self-explanatory. 


8.  PERFORMING  ORGANIZATION  REPORT  NUMBER. 

Enter  all  unique  alphanumeric  report  numbers  assigned 
by  the  performing  organization,  e.g.  BRL-1234; 
AFWL-TR-85-401  7-Vol-21  -PT-2. 

9.  SPONSORING/MONITORING  AGENCY  NAME(S) 

AND  ADDRESS(ES).  Enter  the  name  and  address  of  the 
organization(s)  financially  responsible  for  and 
monitoring  the  work. 

10.  SPONSOR/MONITOR'S  ACRONYM(S).  Enter,  if 
available,  e.g.  BRL,  ARDEC,  NADC. 

11.  SPONSOR/MONITOR  S  REPORT  NUMBER(S). 

Enter  report  number  as  assigned  by  the  sponsoring/ 
monitoring  agency,  if  available,  e.g.  BRL-TR-829;  -215. 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT.  Use 

agency-mandated  availability  statements  to  indicate  the 
public  availability  or  distribution  limitations  of  the 
report.  If  additional  limitations/  restrictions  or  special 
markings  are  indicated,  follow  agency  authorization 
procedures,  e.g.  RD/FRD,  PROPIN,  ITAR,  etc.  Include 
copyright  information. 

13.  SUPPLEMENTARY  NOTES.  Enter  information  not 
included  elsewhere  such  as:  prepared  in  cooperation 
with;  translation  of;  report  supersedes;  old  edition 
number,  etc. 

14.  ABSTRACT.  A  brief  (approximately  200  words) 
factual  summary  of  the  most  significant  information. 

1  5.  SUBJECT  TERMS.  Key  words  or  phrases 
identifying  major  concepts  in  the  report. 

16.  SECURITY  CLASSIFICATION.  Enter  security 
classification  in  accordance  with  security  classification 
regulations,  e.g.  U,  C,  S,  etc.  If  this  form  contains 
classified  information,  stamp  classification  level  on  the 
top  and  bottom  of  this  page. 

17.  LIMITATION  OF  ABSTRACT.  This  block  must  be 
completed  to  assign  a  distribution  limitation  to  the 
abstract.  Enter  UU  (Unclassified  Unlimited)  or  SAR 
(Same  as  Report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited. 


Standard  Form  298  Back  (Rev.  8/98) 


Journal  of  Nuclear  Materials  390-391  (2009)  88-91 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Nuclear  Materials 

journal  homepage:  www.elsevier.com/locate/jnucmat 


Transfer  of  rovibrational  energies  in  hydrogen  plasma-carbon  surface  interactions 

P.S.  Krstic3’*,  E.M.  Hollmannb,  C.O.  Reinhold3,  S.J.  Stuart c,  R.P.  Doernerb,  D.  Nishijimab, 

A.Yu.  Pigarovb 

aOak  Ridge  National  Laboratory,  Physics  Division,  P.O.  Box  2008,  Oak  Ridge,  TN  37831-6372,  USA 
b  University  of  California,  San  Diego,  CA,  USA 

cClemson  University,  Department  of  Chemistry,  Clemson,  SC  29634,  USA 


ARTICLE  INFO 


ABSTRACT 
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We  present  state-of-the-art  molecular  dynamics  (MD)  simulations  of  high-density  plasma-bombardment 
of  carbon,  mimicking  plasma-wall  interactions  at  the  fusion  reactor  first  wall.  Bare  and  hydrogenated 
amorphous  carbon  surfaces  with  temperatures  in  the  range  of  300-800  K  are  bombarded  by  a  distribu¬ 
tion  of  neutral  hydrogen  molecules  representing  well-defined  center-of-mass  and  rovibrational  temper¬ 
atures.  The  MD  simulations  are  benchmarked  against  experiments  in  which  a  heated  carbon  surface  is 
irradiated  with  hydrogen  molecules  from  a  plasma  source.  Comparisons  between  simulations  and  exper¬ 
iment  are  presented  for  the  rovibrational  distributions  upon  reflection  and  the  rotational  and  vibrational 
accommodation  coefficients. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Since  plasma-boundary  physics  encompasses  some  of  the  most 
important  unresolved  issues  in  future  energy  production  for  fusion 
reactors,  there  is  a  strong  interest  in  the  fusion  community  for  bet¬ 
ter  understanding  and  characterization  of  plasma-surface  interac¬ 
tions.  These  interactions  lead  to  surface  erosion  and  particle 
deposition,  which  degrades  fusion  performance  and  produce 
long-term  particle  retention. 

Motivated  by  this  need,  we  have  recently  undertaken  a  series 
of  large-scale  classical  molecular  dynamics  (MD)  simulations  of 
particle-surface  interactions  involving  carbon  surfaces.  Our  previ¬ 
ous  calculations  were  focused  on  modeling  chemical  sputtering, 
implantation,  reflections  and  dissociation  processes  of  hydroge¬ 
nated  carbon  from  well-defined  low-intensity  ion  beam  experi¬ 
ments  [1-3].  In  this  work  we  extend  our  previous  approach  to 
model  a  high-density  plasma-bombardment  environment.  Our 
simulation  results  are  directly  compared  with  carefully  designed 
plasma  irradiation  experiments  performed  at  the  University  of 
California  (San  Diego)  [4—8],  In  Section  2  we  describe  our  MD 
simulations  of  the  plasma-facing  surfaces.  In  Section  3  we  de¬ 
scribe  the  experimental  setup  in  some  detail.  Both  experimental 
and  theoretical  results  and  our  discussion  are  presented  in  Sec¬ 
tion  4. 


*  Corresponding  author. 
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2.  Theory 

Details  of  our  classical  MD  approach,  as  well  as  of  the  bond  or¬ 
der  REBO  potential,  can  be  found  in  previous  publications  (see  [1- 
3]  and  references  therein).  Here  we  describe  only  the  new  aspects 
needed  to  model  plasma  experiments,  i.e.  the  randomization  of  the 
simulated  impact  parameters  to  represent  a  plasma.  Our  previous 
studies  have  shown  [1-3]  that  the  best  agreement  with  experi¬ 
ments  [3]  is  reached  when  the  surfaces  are  prepared  by  bombard¬ 
ment,  closely  mimicking  the  experiments,  until  a  ‘steady-state 
surface’  is  reached.  Applying  this  prescription  to  a  surface  irradi¬ 
ated  by  plasma  requires  that  the  impact  particles  have  the  proper 
distributions  of  energy,  angle,  position,  angular  momentum,  during 
both  the  preparation  of  the  surfaces  and  the  calculation  of  the  par¬ 
ticle  emission  yields. 

2.1.  Sampling  initial  conditions  for  H2  molecules  in  the  plasma 

Clearly,  to  sample  initial  conditions  one  can  separate  the 
internal  degrees  of  freedom  (relative  position  and  momentum 
r,p)  and  the  translational  center-of-mass  (c.m.)  degrees  of  freedom 
(fcm, Pcm)  of  the  diatomic  molecule.  Initial  conditions  can  be  sam¬ 
pled  from  a  probability  density 

p^rem, pcm)  piM(r,p) 

The  phase  space  coordinates  have  to  be  sampled  through  distri¬ 
butions  that  depend  only  on  the  constants  of  motion  of  the  system 
in  order  to  obtain  a  stationary  distribution.  In  our  case  these  are 
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Pcm:  Hint,  L,  where  Hint  is  the  Hamiltonian  of  rovibrational 
motion  and  L  is  the  angular  momentum.  For  the  c.m.  coordinates, 
sampling  is  relatively  simple  and  one  can  use  p^ifcm-Pan)  = 
S(Zcm  -Zo)H(xcm)H(ycm)p^“(pcm)I  where  z0  is  the  initial  z  position 
of  the  c.m.  along  the  surface  normal,  and  ycm  are  sampled  uni¬ 
formly,  and  pan  is  sampled  from  a  shifted  Maxwell-Boltzmann  dis¬ 
tribution  p““(Pcm)  with  mean  velocity  v0  and  translational 
temperature  Tt  and  U(x)  is  a  uniform  distribution  in  x.  For  the 
internal  degrees  of  freedom,  a  stationary  distribution  is  plnt(r,p) 
=  pint(Hint,L)  =  pM(Hint)prot(L)U(L),  where  U(L)  is  a  uniform  distri¬ 
bution  for  the  direction  of  the  angular  momentum  of  the  molecule 
(in  the  unit  sphere),  and  pvib(Hmt),  pmt(L)  are  Boltzmann  distribu¬ 
tions  for  the  ‘vibrational’  and  ‘rotational’  degrees  of  freedom.  We 
first  sample  L  from  prot(L),  use  the  chosen  L  value  into  p„,i,(Hint) 
to  get  (r,pr)  and  then  calculate  p±  =  L/r.  In  order  to  mimic  the 
experiment,  we  use  a  Maxwell  distributions  with  center-of-mass 
temperature  of  Tt  =  800  K,  while  the  internal  temperatures  distri¬ 
butions  were  chosen  with  3000  K  for  rovibrational  motion  (i.e. 
5000  K  for  vibrational  temperatures)  and  1000  K  for  rotational 
motion.  This  lack  of  equilibrium  between  rotational  and  vibra¬ 
tional  temperatures  is  the  best  estimate  of  the  experimental  condi¬ 
tions  in  the  considered  non-equilibrium  plasma  with  short 
confinement  times. 


2.2.  Preparation  of  the  target  surface 

The  target  surfaces  in  our  MD  simulations  are  bare  and  hydro¬ 
genated  amorphous  carbon  (a-C:H).  This  reflects  the  fact  that  the 
graphite  divertor  target  in  the  experiment  was  polycrystalline 
POCA  graphite  surface,  without  a  well-defined  direction  of  the  ba¬ 
sal  plane.  The  experiment  implies  long  plasma  exposure  times  and 
therefore  a  possibility  of  hydrogenation  of  the  surface.  Although 
the  bulk  material  is  known  to  have  a  H:C  ratio  of  0.4  at  saturation 
level,  simulations  have  shown  that  at  low  impact  energies  this  ra¬ 
tio  can  be  substantially  elevated  at  the  interface  [1-3].  Since  the  le¬ 
vel  of  hydrogenation  in  the  experiment  is  not  known,  we  use  both 
bare  (non-hydrogenated)  surfaces  as  well  as  a-C:H  surfaces  ex¬ 
posed  to  varying  levels  of  fluence  with  a  1  eV  beam.  Building 
hydrogenated  surfaces  for  all  possible  plasma  parameters  would 
have  involved  a  bigger  effort.  However,  the  approach  we  undertook 
gave  us  a  clear  answer  on  the  role  of  the  hydrogenation.  In  order  to 
speed  up  the  hydrogenation  the  surfaces  were  prepared  by  cumu¬ 
lative  bombardment  of  an  initial  surface  of  1750  carbon  atoms  and 
700  hydrogen  atoms,  with  carbon  density  of  2  g/cm3  and  at  300  K, 
homogeneously  hydrogenated  at  a  H:C  ratio  of  0.4.  These  were 
then  equilibrated  at  various  temperatures  (300,  400,  600  and 
800  K),  using  a  Langevin  thermostat.  Five  surfaces  were  prepared 
at  fluences  of  100,  200,  400,  800  and  1200  impacts  of  H2  molecules 
on  a  surface  cell  with  dimensions  of  10  A  x  (26.5  A)2,  which  reflect 
various  levels  of  hydrogenations  at  the  surface.  The  two  higher  flu¬ 
ences  are  enough  that  the  surface  has  entered  a  steady-state  re¬ 
gime,  and  surface  composition  changes  little  with  subsequent 
impacts. 


2.3.  Simulation  procedure 

For  each  surface,  defining  a  hydrogenation  level,  and  for  each 
surface  temperature,  the  irradiation  is  repeated  for  1000  random 
independent  trajectories.  The  c.m.,  rotational,  rovibrational  and 
vibrational  energies  and  corresponding  kinetic  temperatures  were 
determined  for  each  reflected  H2  molecule,  and  averaged  over  all 
trajectories.  The  number  of  trajectories  for  each  surface  was  satis¬ 
factorily  large,  resulting  in  less  than  5%  statistical  uncertainty  (one 
standard  deviation)  of  the  calculated  averages. 


3.  Experiment 

The  PISCES-A  reflex-arc  plasma  discharge  [4]  is  used  to  create 
rovibrationally  excited  hydrogen  molecules  which  fill  a  vacuum 
chamber.  These  H2  neutrals  travel  down  a  stainless  steel  side  port 
and  then  enter  a  small  graphite  tube  mounted  in  the  side  port.  An 
electron  beam  is  aimed  down  the  hollow  graphite  tube  to  electron¬ 
ically  excite  a  small  fraction  of  the  hydrogen  molecules  inside  the 
tube.  The  resulting  line  emission  is  then  measured  at  4  positions  in 
the  tube.  The  local  H2  ground  electronic  state  vibrational  and  rota¬ 
tional  temperatures  T and  Twt  at  each  fiber  location  can  then  be 
calculated  with  measured  line  ratios  of  H2  Fulcher  band  emission 
[5,6].  Additionally,  the  H2  kinetic  temperature  Tkin  and  density 
nH2  can  be  obtained  from  line  broadening  and  line  intensity, 
respectively. 

A  detail  of  the  graphite  tube  arrangement  is  shown  in  Fig.  1  (for 
simplicity,  the  adjacent  plasma  discharge  is  not  shown  here).  Col¬ 
limator  lenses  are  used  to  increase  the  amount  of  light  gathered 
into  the  four  fibers.  Heater  wire  wrapped  around  the  graphite  tube 
is  used  for  heating.  The  graphite  tube  temperature  (assumed  to  be 
uniform)  is  measured  with  a  thermocouple  at  the  end  of  the  tube. 
The  graphite  tube  is  made  of  E-294  graphite  from  Electro-Tech 
Machining  with  density  p  =  1 .80  g/cm3  and  average  grain  size  of 
0.10  mm.  The  tube  dimensions  are:  length  L=10cm,  diameter 
D  =  3.8cm,  and  thickness  f  =  0.64  cm.  Prior  to  taking  data,  the 
graphite  was  baked  out  for  4  h  at  1000  °C  to  remove  stored  water. 
The  graphite  was  not  exposed  to  plasma  prior  to  the  experiments, 
so  hydrogen  in  the  walls  arrives  as  neutrals,  not  as  embedded  ions. 

The  electron  beam  is  created  using  a  D  =  1.2  mm  Ta  disk  cath¬ 
ode  emitter.  Typical  beam  energies  are  70  eV;  this  is  chosen  to 
be  well  above  the  threshold  for  H2  electronic  excitation  (~10eV) 
so  that  signal  brightness  is  not  sensitive  to  small  changes  in  beam 
energy.  Beam  currents  (~1  mA)  are  small  and  heating  of  H2  in  the 
tube  due  to  the  beam  is  found  to  be  negligible. 

4.  Results  and  discussion 

Typical  for  all  calculations  performed  here  is  that  the  results  for 
hydrogenated  surfaces  are  different  from  the  results  with  the  bare 
carbon  surface.  Fig.  2  shows  the  average  temperatures  of  various 
modes  of  motion  for  H2  after  reflection  from  the  surface.  The  ini¬ 
tially  similar  values  of  rotational  and  c.m.  temperatures  (1000  K 
and  800  K)  lead  to  very  similar  values  of  the  final  temperatures 
(the  results  for  hydrogenated  surfaces  lie  within  the  same  depicted 
shaded  area).  However,  the  reflected  H2  is  typically  ‘hotter’  both 
rotationally  and  translationally  for  the  hotter  surfaces,  and  hotter 
for  the  bare  than  for  the  hydrogenated  surfaces,  becoming  cooler 
with  increased  hydrogenation  (arrow  in  Fig.  2).  On  the  other  hand, 
the  vibrational  temperature  is  significantly  decreased  in  the  colli¬ 
sion,  though  no  clear  trend  is  found  with  the  hydrogenation  level. 
These  observations  imply  that  the  reflection  at  the  surface  is 
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Fig.  1.  Close-up  of  graphite  tube  mounted  in  side  port. 
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Fig.  2.  Post-collisional  temperatures  of  reflected  H2  as  a  function  of  the  surface 
temperature.  The  curves  for  hydrogenated  surfaces  are  shown  by  shaded  areas 
between  dashed  lines,  while  the  bare  surface  results  are  displayed  by  the  thick  solid 
lines  with  open  symbols.  The  center-of-mass  translational  temperature  is  given  by 
the  short  dashed  line  with  solid  squares.  The  thin  doted  lines  represent  temper¬ 
atures  of  the  impacting  H2  molecules.  The  arrows  indicate  the  direction  in  which 
the  results  change  with  increasing  level  of  hydrogenation. 

accompanied  by  the  conversion  of  part  of  the  vibrational  (and 
surface)  energy  into  translational  and  rotational  energy.  The  mole¬ 
cule  is  also  cooled  down  rovibrationally  in  the  collision,  and  this 
effect  increases  with  hydrogenation,  similarly  to  the  rotational 
temperature. 

Fig.  3  shows  the  change  in  energy  of  the  various  modes  of  H2  in 
the  collision  with  the  surface,  consistent  with  the  discussion  in 
Fig.  2.  The  positive  sign  of  the  energy  change  for  the  rotational 
and  c.m.  energy  is  consistent  with  a  negative  accommodation  coef¬ 
ficient  for  these  modes  of  motion  in  the  simulated  system,  which  is 
not  confirmed  by  our  experiment,  as  illustrated  in  Fig.  5.  Quite  pos¬ 
sible  is  that  the  coupling  (mixing)  of  vibrational  and  rotational 
modes  of  motion,  obtained  here  classically,  is  quantum-mechani- 
cally  suppressed  due  to  a  large  energy  difference  of  the  vibrational 
energy  levels  and  the  low  level  of  vibrational  excitation.  Therefore, 
our  simulation  results  should  be  accepted  with  reserve  when  con¬ 
sidering  separate  rotational  and  vibrational  modes,  and  considered 
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Fig.  3.  Energy  transfer  of  the  different  modes  of  motion  from  H2  molecules  in  a 
collision  with  bare  (solid  thick  line  with  symbols)  and  hydrogenated  surfaces  of 
various  temperatures. 


Fig.  4.  The  rovibrational  accommodation  coefficient  for  bare  and  hydrogenated 
surfaces,  at  various  surface  temperatures.  The  description  of  the  lines  is  the  same  as 
in  Fig  2. 

mainly  for  the  conclusions  on  their  sum,  the  rovibrational  energy, 
which  is  here,  unlike  separate  rotational  and  vibrational  energies,  a 
conserved  quantity.  The  accommodation  coefficient  arovib  = 
(Tf  —  Ti)/(TS  —  Tt),  for  rovibrational  motion  upon  collision  of  H2 
with  the  surface  is  shown  in  Fig.  4.  Tf  i  are  final  and  initial  rovibra¬ 
tional  temperatures  of  the  molecule,  while  Ts  is  the  surface  temper¬ 
ature.  It  shows  a  trend  of  increase  with  the  level  of  hydrogenation, 
consistent  with  higher  energy  transfers  in  collisions  with  lighter 
particles.  Due  to  the  large  difference  between  the  vibrational  and 
rotational  temperatures,  arovib  is  dominated  by  the  largest  temper¬ 
ature  (vibrational).  Interestingly,  calculated  accommodation  coeffi¬ 
cients  are  within  the  uncertainty  margins  of  the  experimental 
results  below. 

To  interpret  the  measurements,  Monte  Carlo  modeling  of  H2 
neutral  trajectories  is  used  [7],  The  free  parameters  in  the  modeling 
are  the  temperatures  Tvib ,0  and  Trol0  at  the  tube  entrance  and  the 
accommodation  coefficients  arot  and  ocvib.  Wall  collisions  are  as¬ 
sumed  to  be  diffuse.  The  kinetic  temperature  is  assumed  to  be  al¬ 
ways  equal  to  the  rotational  temperature,  Tkin  —  Trot t  as  is  seen  in 
the  experiments.  Because  Tvib  is  found  to  be  much  larger  (factor 
of  5-10)  than  Trot,  direct  mixing  between  Trot  and  T„ih  in  wall  colli¬ 
sions  is  assumed  to  be  small  and  is  neglected.  Neutral-neutral  col¬ 
lisions  and  rovibrational  energy  transfer  between  vibrational 
energy  and  kinetic/rotational  energy  in  these  collisions  is  included 
in  the  Monte-Carlo  modeling.  Neutral  density  is  obtained  from  the 
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Fig.  5.  Monte-Carlo  fits  (lines)  of  experimental  data  (symbols)  for  (a)  Trot  vs  axial 
position  and  (b)  Tvib  vs  axial  position  (see  text  for  a  description  of  the  lines). 
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Fig.  6.  Measured  (a)  rotational  accommodation  coefficients  and  (b)  vibrational 
accommodation  coefficients  as  a  function  of  surface  temperature  Ts.  Dashed  lines 
are  fits  through  the  data.  Shaded  symbols  are  for  H2  and  white  symbols  are  for  D2. 

brightness  of  the  H2  lines;  this  is  found  to  depend  dominantly  on 
the  temperature  of  the  graphite  tube.  The  standard  definition  of 
the  accommodation  coefficient  is  used  [8],  i.e.  after  a  wall  collision 
the  new  rotational  temperature  is  Trot  =  Tro,  -  ocrot(Twt  -  Ts),  where 
Ts  is  the  surface  temperature.  An  analogous  definition  is  used  for 
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Fig.  5  gives  an  example  of  Monte-Carlo  fits  to  the  measured  Tmt 
and  Tvib.  The  experiment  shown  here  was  for  room-temperature 
graphite.  The  free  parameters  Troti0,  Tvib 0,  ccwt,  and  were  varied 
to  obtain  a  best  fit  shown  by  the  solid  lines.  The  error  in  the  accom¬ 
modation  coefficients  was  estimated  by  finding  the  highest  and 
lowest  values  of  arot,  and  txvib  which  would  still  fit  the  data  within 
the  error  bars;  these  extreme  fits  are  shown  by  the  dashed  curves 
in  Fig.  5.  Error  bars  on  Trot  and  Tvib  data  are  estimated  from  the 


quality  of  the  single-temperature  Boltzmann  fits  to  the  line 
populations. 

Fig.  6  shows  the  values  of  the  accommodation  coefficients  a.m, 
and  avib  obtained  from  various  experimental  runs  (e.g.  each  point 
in  Fig.  6  corresponds  to  a  fit  such  as  shown  in  Fig.  5).  Overall,  the 
rotational  accommodation  aw,  can  be  seen  to  be  somewhat 
2  -  lOx  larger  than  the  vibrational  accommodation  a.vib.  Within 
the  scatter  of  the  data,  D2  accommodation  (white  symbols)  is  sim¬ 
ilar  to  H2  accommodation  (shaded  symbols).  A  slight  increase  in 
accommodation  with  increasing  surface  temperature  is  suggested 
by  the  data. 
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